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increases for bigger velocity fluctuation amplitude and is associated to anisotropy of the studied 
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PACS numbers: 05.20.Jj, 47.27.Jv 



I. INTRODUCTION 

The nonextensive statistics developed by Tsallis Q] has 
motivated a phenomenological approach 0- 0- 01 which 
was recently used |5|, |6( to describe Lagrangian accelera- 
tion of fluid particle in developed three-dimensional tur- 
bulent flow within the framework of Langevin type equa- 
tion; see also @,H,E1E|. Some recent stochastic particle 
models and refinements of this technique 0,0, El were 
reviewed in Ref. [T^| . The one-dimensional Langevin 
toy models of Lagrangian turbulence lacks of phys- 
ical interpretation, e.g., of short term dynamics, or 
small-scale and large-scale contributions, in the context 
of three-dimensional incompressible Navier-Stokes equa- 
tion, which is generally assumed to be capable to de- 
scribe fully developed turbulence. A deep physical anal- 
ysis of this problem has been made by Gotoh and Kraich- 
nan fl5^ . 

Recently [H EH we have shown that the one- 
dimensional Laval-Dubrulle-Nazarenko (LDN) type toy 
model [T7L ITsI ] of the acceleration dynamics with the 
model turbulent viscosity v t and coupled delta-correlated 
Gaussian multiplicative and additive noises is in a good 
agreement with the recent hi gh-p recision experimental 
data on acceleration statistics |2(J, |21j . Particularly, 
we have demonstrated that the predicted contribution 
to fourth order moment of acceleration does peak at the 
same values as the experimental curve, in contrast to 
predictions of most of the other stochastic particle mod- 
els 0, 0, ^3 ^M- We note however that characteristic 
parts of the acceleration distribution are the core and 
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tails, while the intermediate range of accelerations which 
is responsible for positions of the peaks may be not ro- 
bust. 

Below we discuss at some length on foundations of the 
present model, and focus on modeling acceleration statis- 
tics conditional on velocity fluctuations. 

The original one-dimensional LDN model was formu- 
lated both in the Lagrangian and Eulerian frameworks for 
small-scale velocity increments. It is based on a stochas- 
tic kind of Batchelor-Proudman rapid distortion theory 
(RDT) approach to the three-dimensional Navier-Stokes 
equation and Gabor transformation • Small scales are 
separated and assumed to be stochastically distorted by 
much larger scales. Such nonlocal interactions which can 
be viewed as elongated triads are featured by the three- 
dimensional LDN model, to which the one-dimensional 
model makes an approximation, particularly via mimick- 
ing relationship between stretching and vorticity. Local 
interactions are modeled by the turbulent viscosity which 
is taken due to the renormalization group approach. 

In the present paper, we use the Lagrangian formula- 
tion of the one-dimensional LDN model, which is char- 
acterized by a simple structure, as an ansatz to formu- 
late Langevin type equation for the component of fluid 
particle acceleration in statistically homogeneous and 
isotropic developed turbulence. 

In general the Langevin type equation which we em- 
ploy here contains wellknown terms, namely, time deriva- 
tive of the variable, an additive noise, and a nonlinear 
drift term with a multiplicative noise; see Eq. below. 
The noises represent large scales and are treated indepen- 
dent on small scales. In the Lagrangian frame, they are 
taken in the simplest realization, Gaussian white-in-time 
with zero mean, and are assumed to be delta-correlated 
to each other. Such choice of the noises correspond to 
the zero-time approximation of small time-scale corre- 
lations viewed along the Lagrangian trajectory of fluid 
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particle. Direct account for finite time correlations in a 
Langevin type equation is important but this is beyond 
the scope of the present study, in which we concentrate 
on the longtime behavior and qualitative analysis. Long- 
time dynamics of large scales is ignored in this particular 
case of the LDN model. 

In the conventional approach, wandering of the tracer 
particle can be described within the framework of gener- 
alized Brownian like motion, with the acceleration cho- 
sen to be the random variable depending on time t (a 
stochastic process) . The statistical isotropy allows one to 
decompose the acceleration into the 'transverse' and 'lon- 
gitudinal' components with respect to the corresponding 
velocity that reduces formally the analysis to effective 
two-dimensional consideration of the essentially three- 
dimensional system. The isotropy is assumed to be lo- 
cal, in the neighbor of fluid particle. In one-dimensional 
models one considers one of the two components not cou- 
pled to each other. Joint dynamics of the components of 
acceleration is of much interest and can be considered 
elsewhere. 

It should be emphasized here that the usual Eulerian 
framework (fixed probe) represents the structural view- 
point which is extensively studied in the literature while 
the Lagrangian framework is associated to stochastic dy- 
namical consideration, a statistical viewpoint. These 
two basic approaches are different both in their theoret- 
ical formulations and experimental technique, and com- 
pliment each other. In the Eulerian frame, the accel- 
eration, a, = dvi/dt = dfVi + VkdkVi, is expressed in 
terms of the velocity field and its temporal and spa- 
tial derivatives, while in the Lagrangian frame the ac- 
celeration is ai = dtVi, in terms of the Lagrangian ve- 
locity, Vi = dtXi, of fluid particle with the coordinates 
Xi(xk(0),t); i,k = 1,2,3. In the present paper we will 
use the Lagrangian framework. 

Measurement of the time series Xi(t) for individual 
tracer particle in a turbulent flow to a high resolution re- 
quires very high speed imaging sensors, and allows one to 
compute its instantaneous Lagrangian velocity and La- 
grangian acceleration |l9j . Experimental data are col- 
lected for such individual particles moving in the ap- 
proximately isotropic and homogeneous turbulent flow 
domain. 

With delta-correlated noises the Langevin type model 
of one of the two components of acceleration is Markovian 
(no memory effects) so that the well established Fokker- 
Planck approximation can be used to link the dynamical 
framework to a statistical (PDF) approach. Thus, one 
can study stationary solution of the associated Fokker- 
Planck equation for the acceleration one-point probabil- 
ity density function P(a, t) under the assumption of bal- 
ance between the energy injected by driving forces at 
large scales and the energy dissipated by viscous pro- 
cesses at small scales (a statistically steady state). 

In the case when stationary probability distribution 
can be found exactly one can make further analysis with- 
out referring to dynamics. The resulting stationary prob- 



ability density function of the component of acceleration, 
P(a), contains some free parameters inherited from the 
Langevin type equation. This PDF is the main predic- 
tion of the model which can be directly fitted to experi- 
mental data or direct numerical simulation (DNS) of the 
Navier-Stokes equation. 

The Random Intensity of Noise (RIN) approach [Til 
[T^ | makes a simple extension of the above Langevin type 
model of the Lagrangian acceleration. The main idea of 
this approach is to use the presence of two well sepa- 
rated characteristic time scales (Kolmogorov time scale 
and the integral time scale) of the system and assume 
that parameters entering the resulting acceleration PDF, 
such as intensities of the noises, are not constant but 
fluctuate at large time scale and depend on Lagrangian 
velocity fluctuations. 

Effectively, we approximate the time evolution of ac- 
celeration by separating fast and slow time varying parts, 
and thus account for large time scales which are usually 
ignored. This improves the delta-correlation approxima- 
tion of the noises adopted in the one-dimensional LDN 
type model mentioned above, provided that character- 
istic longtime fluctuations along the particle trajectory 
occur mainly at time scales much bigger than the dissi- 
pative time scale, typically at the integral time scale. 

In general, one can assume hierarchy of a number of 
characteristic time scales associated to the discrete cas- 
cade picture with characteristic times of the flow modes 
of descending eddies. However, in the present paper we 
simplify the consideration by extending the large time 
scale up to the Lagrangian integral time scale in order to 
make it more analytically tractable, in accord to the pres- 
ence of two basic characteristic scales in the Kolmogorov 
1941 (K41) picture of fully developed turbulence, and the 
recent Lagrangian experimental data p2^ | . 

In the context of stochastic equation the delta- 
correlated multiplicative noise provides short-time accel- 
eration bursts the origin of which is thought to be due to 
the presence of very intense vortical structures relatively 
slowly varying in time. These structures make an essen- 
tial contribution to the acceleration statistics, namely, 
to the tails of acceleration probability density function, 
making it highly non-Gaussian. Accounting for longtime 
fluctuations of the parameters corresponds to an account- 
ing for fluctuations of intensity of the vortical structures, 
intensity of noisy incoherent background, and their cross 
correlation. The key point is then to identify distribu- 
tions of the parameters. We will discuss it below and 
further in Sec. IIIII 

It is worthwhile to mention that accounting for the 
fluctuating parameter (3 in the chi-square Langevin model 
characterized by the sole delta-correlated Gaussian ad- 
ditive noise |(j was recently found [14| to yield a sta- 
tionary probability density function of the same power 
law form as in simple Langevin model characterized by 
the delta-correlated Gaussian additive and multiplicative 
noises with constant parameters and linear drift term. In 
other words, the chi-square distributed fluctuations of j3 



3 



appears to be mimicking the presence of the multiplica- 
tive noise in this particular case. 

In the first approximation the velocity statistics in 
the Lagrangian frame is taken as usually stationary and 
Gaussian, partially because one can easily deal with it an- 
alytically. Due to the recent Lagrangian experiments 0] 
the distribution of Lagrangian velocity is approximately 
Gaussian for both the x and z components (the flatness 
is 3.2 and 2.8 respectively as compared with 3 for a Gaus- 
sian) and their characteristic time variation is of the or- 
der of Kolmogorov time, yet velocity fluctuations exceed 
root-mean square (rms) velocity that is usually associ- 
ated to large scales. This can be understood as a mani- 
festation of the importance of nonlocal (inter-scale) inter- 
actions when considering small scales. It is important to 
note here that the Lagrangian velocity autocorrelation, 
as well as the autocorrelation of absolute value of velocity 
increments in time, was found to cross zero at the inte- 
gral time scale, while the full signed velocity increments 
in time decay at the Kolmogorov time scale j2^| . 

Under the assumption that the parameters are given 
independent random variables (large scales are weakly 
affected by small scales in a three-dimensional high- 
Reynolds-number flow since the former are local in the 
wave number space) characterized by stationary statis- 
tics, the stationary probability density function obtained 
from the Fokker-Planck equation makes a sense and is 
treated as the distribution conditional on small scale ve- 
locity fluctuations through assumed dependencies of the 
parameters on u. 

By this way the RIN approach enables one to study ac- 
celeration statistics conditional on velocity fluctuations, 
P(a\u) — P(a|Parameters(u)). The velocity fluctuations 
(large scales) are assumed to be decoupled from the accel- 
eration (small scales) at high Taylor microscale Reynolds 
numbers, R\ > 500 |19| . 

It should be noted that such an approach is in agree- 
ment with the Heisenberg-Yaglom picture of developed 
turbulence which relates statistics of the fluid particle 
acceleration to statistics of velocity fluctuations on the 
basis of K41 scaling theory, with pressure gradient con- 
tribution to the acceleration variance strongly dominat- 
ing over that of viscous forces, in the inertial range. The 
long-standing Heisenberg-Yaglom scaling, (a 2 ) ~ u 9 / 2 , 
where u is rms velocity, was recently confirmed to a high 
accuracy for seven orders of magnitude in (a 2 ), and was 
found to be broken for R\ < 500 pjj due to increasing 
coupling of the acceleration to large scales of the flow. 

Also, it should be mentioned that a similar approach, 
with the variance of intermittent variable viewed phe- 
nomenologically as a parameter which follows log-normal 
distribution, was considered by Castaing, Gagne, and 
Hopfinger [2jj . From this point of view RIN models can 
be referred to as Castaing type models. However, the 
Castaing model does not refer to a stochastic dynamical 
approach, which is an important ingredient of the La- 
grangian modeling, and does not relate the fluctuating 
parameter to velocity fluctuations. 



In the RIN approach the marginal (i.e., unconditional) 
stationary probability density function P(a) is found sim- 
ply by integrating out u in the conditional distribution 
P(a\u) with (Gaussian) distributed u. This procedure 
requires prior determination of dependencies of the pa- 
rameters on velocity fluctuations u. 

Particularly, in the RIN approach to the LDN type 
model the assumption that the additive noise intensity a 
depends on absolute value of velocity fluctuations u in an 
exponential way, a ~ e'"', was found [l4j to imply a set 
of conditional probability density functions P(a\u) and 
the conditional acceleration variance (a 2 \u) which are in 
a good qualitative agreement with the recent experimen- 
tal data on the conditional statistics of the x component 
of acceleration obtained by Mordant, Crawford, and Bo- 
denschatz for the normalized velocity fluctuations 
^/(it 2 ) 1 / 2 ranging from to 3. This issue will be dis- 
cussed further in Sec. IIIII 

The effect of nonzero cross correlation A between the 
additive and multiplicative noises (which models the rela- 
tionship between small-scale stretching and vorticity [l7| 
and can be seen as a small skewness of the probability 
density function of the 'longitudinal', i.e., pointed along 
the Lagrangian velocity, component of acceleration) has 
been studied in Ref. [Ig for an illustrative purpose. The 
experimental unconditional acceleration probability den- 
sity function reveals small skewness that has been fitted 
by using the value A = —0.005. 

The observed very small skewness of the probability 
density functions for the x component of acceleration 
can be assigned to anisotropy of the studied R\ = 690 
flow [l6l l2lj rather than to the effect of correlation be- 
tween stretching and vorticity. Particularly, it was found 
that, in addition to the dependence a — a(u), the param- 
eter A should also depend on velocity fluctuations to meet 
the experimentally observed appreciable increase of the 
conditional mean acceleration, (a\u), with the increase 
of |u| l2ll |. The form of functional dependence \(u) was 
assumed to produce a weak effect but it was not specified. 

In the present paper we fill this gap and demonstrate 
that (rather strong) exponential dependence, A ~ e'"', is 
in a good qualitative agreement with the experimental 
data. This fit is of a pure illustrative character since the 
observed nonzero conditional mean acceleration is due 
to the anisotropy of the studied flow. We perform this 
fit to study influence of variation of the cross correla- 
tion parameter on the shape of conditional acceleration 
distribution and conditional acceleration variance. 

The layout of the paper is as follows. In Sec. [H] we 
outline results of the one-dimensional LDN type model of 
the component of acceleration. In Sec. IIIII we study the 
conditional mean acceleration using the RIN extension 
of the LDN type model, with the additive noise intensity 
a ~ e' u ', the cross correlation parameter A ~ e' 11 ', and 
the other parameters of the model fixed. We summarize 
the obtained results in Sec. IIVI 
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II. THE LDN TYPE LANGEVIN MODEL 

In this Section, we give only a brief sketch of the LDN 
model and refer the reader to Refs. [Til Il7| for more 
details; see also recent paper [24| for the stochastic RDT 
approach. 

We use the exact result for probability density func- 
tion of the LDN type model obtained as a stationary 
solution of the Fokker-Planck equation associated to the 
one-dimensional Langevin equation for the component of 
acceleration 0,0], 



dta = (£ - v t k 2 )a + a±. 



(1) 



This equation is a Lagrangian description in the scale 
space, in the reference frame comoving with the wave 
number packet. This toy model can also be viewed 
as a passive scalar in a compressible one-dimensional 
flow yj|. Here, £ and a± model stochastic forces in the 
Lagrangian frame and are chosen to be Gaussian white- 
in-time noises, 

«(*)) = o, (t(tW)) = 2D6(t-n 

(a ± (t)) = 0, (a ± (t)a x (t')) = 2aS(t - t'), (2) 

where the averaging is over ensemble realizations. The 
probability density function was calculated exactly 0] 
and appeared to be of rather complicated form despite 
the simplicity of the starting equation Q)-©- It is given 
by 



P(a) 



(3) 



Cexp[-v t k 2 /D + F(c) + F(-c)] 
{Da 2 ~2Xa+a) 1 / 2 {2Bka + v t k 2 ) 2BXk / D2 

where we have denoted 

cifc 2 2D 3 
^ ~ 2c 2 D 2 c 11 C!C 2 {c~ Da + X) 
x (B 2 (X 2 + cA - Da)a + c{Dv 2 k 2 + c 2 i/ t ))], (4) 

c = -iy/Da- A 2 , v t = ^v 2 + B 2 a 2 /k 2 , (5) 

ci = B 2 (AX 3 +AcX 2 -3Da\-cDa) + D 2 (c+X)i'^k 2 , (6) 

c 2 = ^B 2 (2X 2 + 2cX - Da)k 2 + D 2 v 2 k A , (7) 

and C is normalization constant. 

Without loss of generality one can put, in a numerical 
study, k — 1 and the additive noise intensity parame- 
ter a — 1 by appropriate rescaling of the multiplicative 
noise intensity D > 0, the turbulent viscosity parameter 
B > 0, the kinematic viscosity Vq > 0, and the cross cor- 
relation parameter A 0], and make a fit of P(a) to the 
experimental data. 

The particular cases B = and v$ = with A put to 
zero were studied in Ref. [l4|. As one can see from Q, 
the real parameter A is responsible for an asymmetry of 
the distribution with respect to a — > —a. For A = the 
distribution is symmetrical. Note that Da — X 2 makes 



a special case, and P(a) was found to be well defined 
for \a\/{a 2 ) 1 / 2 < 60 in the practically interesting case of 
negative A with |A| <SC a and |A| <C D. One observes also 
a rather nontrivial dependence of P(a) on a, X and other 
parameters through F(c) given by Eq. (@J. 



III. THE CONDITIONAL MEAN 
ACCELERATION 

In this Section, we use the exact distribution as a 
starting point. 

Since the experimental unconditional distribution P(a) 
and the experimental conditional distribution P(a\u) at 
u = are a ppr oximately of the same stretched expo- 
nential form |2l| we use the result of our fit 0] of the 
probability density function @ to the unconditional dis- 
tribution P(a) [20( | . This implied the following set of 
values of the real parameters: 



k= 1, a = l, D 
v Q = 2.910, A = 



1.100, B = 0.155, 
-0.005, C = 3.230. (8) 



It should be stressed that a fit to the experimental 
conditional distribution P{a\u) at u = would yield a 
different particular set of values of the parameters. The 
above fit is however justified as a first step, as we are 
mainly interested in a qualitative analysis. Also, the rea- 
son that we do not use a fit to the experimental P(a\u) 
at u = presented in Ref. [2lJ is an illustrative character 
of the curve and that the shown range, \a\/{a 2 ) x / 2 < 15, 
is too small to capture information encoded in the long 
tails that is essential for an accurate determination of the 
fit parameters. 

Following the RIN approach, we will assume that the 
parameters a and A entering © are stochastic and de- 
pend on velocity fluctuations u. In the present paper, 
the remaining parameters, k, -D, B, and isq, are taken to 
be fixed at the fitted values given in (JHJ). 

The exponential form of a{u) has been studied in 
Ref. ^(| and was found to be relevant from both the 
theoretical and experimental points of view. Namely, for 
Gaussian distributed u the positive parameter a is log- 
normally distributed variable that corresponds to Kol- 
mogorov 1962 (K62) refined theory which assumes log- 
normal distribution of the stochastic energy dissipation 
rate per unit mass, e, to which we relate a due to the 
K62 universality hypothesis, in the inertial range. 

The following remark is in order. More precisely, in- 
stead of the stochastic energy dissipation rate it seems 
reasonable to use here the stochastic energy flux through 
the surface of the domain of a given spatial scale. The 
fluctuating energy flux is related to the nonlinear term in 
the Navier-Stokes equation which dominates in the iner- 
tial range for high-Reynolds-number flow, while the fluc- 
tuating energy dissipation rate is associated specifically 
to the dissipative scale. Under the stationarity condition 
(the mean energy dissipation e is equal to the injected 
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energy), the time variation of the energy contained in 
the given scale is defined by the fluctuating energy flux 
minus fluctuating part of the energy dissipation rate at 
the given scale. 

From the experimental point of view, the exponen- 
tial form of a(u) leading to the log-normal R1N model 
yields the acceleration probability density function with 
one free parameter which is in a good agreement with 
tails of the experimental acceleration distribution 0, . 

Clearly, this approach builds first approximation since 
the stochastic energy dissipation rate is known to fol- 
low log-normal distribution only approximately. Also, we 
note that the stochastic energy dissipation rate is posi- 
tive while velocity fluctuations u are not. Hence one is 
motivated to seek for specific model forms of the depen- 
dence between two stochastic variables (monotonic Borel 
function), among which we choose an exponential form 
for the function a(u) that means the same form for the 
function e(u) as we will see below. 

Statistical properties of the stochastic energy dissipa- 
tion rate and of the velocity fluctuations are not the same. 
However, they are related to each other, particularly due 
to the wellknown K62 similarity hypothesis. We assume 
the relationship In a ~ In e in a statistical sense which 
is evidently insensitive to the details related to a power 
law functional dependence of a on e, i.e.. insensitive to 
this particular type of nonlinearity [5J, 1 1 3| . This can be 
viewed a manifestation of the universality. The choice 
of an exponential dependence for a(u) and a power law 
dependence for a(e) implies the relationship e ~ e u , by 
which we model the nonlinear relationship between e and 
u. It should be stresses that only absolute value of u 
contributes the acceleration probability density function, 
for normally distributed u. The framework for dealing 
with more general situation has been recently proposed 
inRef. 0. 

From the phenomenological point of view, the expo- 
nential form, a(u) ~ e'"', was found to provide appre- 
ciable increase of the conditional acceleration variance 
(a \u) with increasing |u| that meets the experimental 
data 

Guided by the above observations the simplest choice 
is to try an exponential dependence for X(u), in a phe- 
nomenological way. Particularly, in the present paper we 
take 

a(u) = e W, A(u) = -0.005e 3|u| , (9) 

which recover the values given in Eq. {HJ at |u| =0. The 
procedure is to refine the guess on X(u) depending on the 
result. 

Of course, the determination of functional forms for 
both a(u) and X(u) is essentially (K62) phenomenological 
but in general this approach is justified from the turbu- 
lence dynamics and allows us to deal with the conditional 
statistics of the fluid particle acceleration, which exhibits 
a quite nontrivial behavior. It should be emphasized that 
the very model is justified by the Navier-Stokes equation 
based LDN model and the very dependence of the 
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FIG. 1: The conditional acceleration probability density 
function P(a\a(u), A(it)) given by for k = 1, a = e' u ', 
D = 1.100, B = 0.155, uq = 2.910, A = -0.005e 31111 ; 
u = 0,0.25,0.50,0.75,1.00,1.19 in the rms units. The top 
curve is conditional on u — while the bottom curve is con- 
ditional on u = 1.19. x — a/(a 2 ) 1 ^ 2 denotes normalized ac- 
celeration. 



additive noise intensity a (and the cross correlation A) 
on u is specific to the LDN model, in which the addi- 
tive noise was found to depend on small-scale velocities 
coupled to large-scale velocities. 

The normalization constant C in will be calculated 
for each value of u, while C — 3.230 corresponds to the 
case u = 0. The conditional probability density function 
is thus given by 10 treated in the form P(a\a(u), A(it)) 
with the parameters defined by JSJ and ©. 

We are now in a position to investigate numerically 
how the variation of u affects the conditional statistical 
properties of the component of acceleration of fluid par- 
ticle. 
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FIG. 2: The conditional mean acceleration as a func- 
tion of standardized velocity fluctuations u. Triangles: 
<a|u)/(a 2 |0) 1/2 for a = e |u| and A = -0.005e 3|u| ; k = 1, 
D = 1.100, B = 0.155,^ = 2.910. Boxes: the experimental 
data on (a\u) / (a 2 ) 1 ^ 2 |21|. it is given in rms units. 



For the normalized velocity fluctuations values 
u/(u 2 } 1/2 = 0, 0.25, 0.50, 0.75, 1.00, 1.19 the dis- 
tributions P(a\a(u) 1 X(u)) and mean accelerations 
(a|u)/(a 2 |0) 1 / 2 are shown in Figs. ^andE! For conve- 
nience the curves in Fig. ^ were shifted by using repeated 
factor 0.1. 

One observes a good qualitative correspondence of 
the obtained conditional distributions plotted in Fig. 2] 
with the experimental curves (Fig. 6a in [2l|])- Note 
that both the variance and skewness of the distribu- 
tion P(a\u) — P(a\a(u), X(u)) increase for bigger velocity 
fluctuations While the increase of the variance [re- 
lated to the increase of ot{u)\ is readily seen, the increase 
of the skewness [related to the increase of X(u)} results 
in rather small change of the shape of distribution de- 
spite the fact that |A| increases by about two orders of 
magnitude, from |A| =0.005 to 0.18. This change could 
be however readily seen in a plot of the contribution to 
fourth order conditional moment, a 4 P(a\u), as it is, e.g., 
illustrated in Fig. 4 of Ref. [jfjj. Notice that the tails 
of the predicted conditional acceleration distributions in 
Fig. ^ remain almost the same while the central peaks 
become weaker for bigger amplitudes of velocity u. 

The obtained conditional mean acceleration 
(a\u)/(a 2 \0) 1/2 plotted in Fig. H is also m a good 
qualitative agreement with the experimental dependence 
(a\u)/(a 2 y/ 2 (Fig. 6b in fift). The mean acceleration 
is evidently zero for symmetrical distribution (A = 0). 
One observes an appreciable increase of the mean 
acceleration for bigger velocity fluctuations We note 
that the experimental dependence of (a\u) / {a 2 ) 1 / 2 on 
u exhibits some asymmetry with respect to u — ► —u. 
This feature is not captured by the model since we have 
chosen A ~ e' u which is a symmetric function. 

It should be stressed that the obtained fit for the mean 
acceleration is of purely illustrative character. It demon- 
strates effects produced by nonzero parameter A which 



depends on velocity fluctuations u, on a qualitative level. 
Our primary purpose was to quantify to which extent 
variation of A affects the conditional variance of accel- 
eration. This effect has been found ignorably small. 
Statistically homogeneous and isotropic turbulence is 
characterized by zero mean acceleration for any compo- 
nent in the laboratory frame of reference. The observed 
nonzero conditional mean acceleration was claimed to re- 
flect anisotropy of the studied von Karman flow [2lJ ■ It is 
however interesting to note that DNS also reveals slight 
departures from zero. Probably this is due to approxi- 
mate character of the isotropy used in the DNS. 

The reason that we make this illustration in the present 
paper is that for the components of acceleration which are 
aligned to trajectory of fluid particle, some skewness of 
the acceleration distributions may be present. Namely, 
for the component pointed transverse to the velocity, A 
is zero by construction while for the component pointed 
along the velocity it is nonzero. We remind that the ac- 
celeration is calculated due to Lagrangian longitudinal 
velocity increments, u(t + r) — u(t) = ra(t), in the dissi- 
pative time-scale r. It is known that Lagrangian longitu- 
dinal velocity structure functions of odd order are small 
but not zero for inertial time-scales so that small skew- 
ness of the velocity increments PDF should be observed. 
We expect that this skewness persists in the dissipative 
range of time-scales that implies (small) skewness of the 
acceleration distribution for the corresponding compo- 
nent. 

Below we make some remarks regarding numerics. 

(i) We have restricted the present numerical study 
by the upper value ^/(m 2 ) 1 / 2 = 1.19 because the dis- 
tribution P(a\a(u), X(u j) given by turns out to be 
ill-defined (discontinuous drop appears at some positive 
value of a) for bigger normalized velocity fluctuation val- 
ues. If this is not due to a failure of the used numerical 
procedure at big |A| in intermediate calculations, it may 
mean that some adjustment of the parameters k, D, B, 
or vq is required to get well-defined P(a\a(u), X(u)). Note 
that it is not presumably well-defined in the entire do- 
main of allowed parameters values due to the presence 
of logarithm in J3J and imaginary terms. Particularly, 
we have found that certain increase of the value of B 
or vq implies well-defined conditional distributions for all 
\u\j (u 2 ) 1 / 2 up to 1.5 but again these become not contin- 
uous for bigger values. Formally, the appearance of dis- 
continuous drop of P(a\u) at some positive a for big |u| 
leads to a steeper increase (saturation) of the conditional 
mean acceleration {a\u} for big |u| which is however not 
observed up to [^[/(m 2 ) 1 / 2 = 2.5, and the experimental 
conditional distributions do not exhibit such a behavior 
up to \u\/ '(u 2 ) 1 / 2 — 3. In the physical context, the condi- 
tion that one should avoid the emerging ill-defined char- 
acter can be understood as that in addition to a(u) and 
X(u) some of the other parameters should depend on u in 
certain way to provide well-defined distributions for any 
u. This is important in various aspects, e.g., to provide 
the integration over u £ [— oo, +oo] to get the marginal 
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distribution P(a). This issue is of much importance to 
the present formalism but it is beyond the scope of the 
present paper and can be considered elsewhere. Another 
possible reason of the ill-defined character is that the con- 
ditional mean acceleration should be very small or zero 
that requires much steeper increase of the parameter A 
for bigger \u\. 

(ii) In Fig.[21we have used (a 2 )!)) 1 / 2 instead of (a 2 ) 1 / 2 
for normalization of the conditional mean acceleration 
(a\u). This makes a change in the overall constant fac- 
tor which is obviously not of much importance in the 
qualitative analysis made in the present paper. For ex- 
ample, smaller value of (a^O) 1 / 2 would shift the whole 
curve (triangles) in Fig. up in the vertical direc- 
tion. A direct numerical calculation of the mean square, 
( fl2 ) = I-oo da a2 f-oo du P{ a W{u), \{u))g(u), where g(u) 
is Gaussian distribution, with the above set up reveals 
slow divergency which is associated to the ill-defined 
character of P(a\a(u), A(u)) at big \u\ mentioned above. 

IV. SUMMARY 

(i) Using the RIN approach to the LDN type one- 
dimensional Langevin model of fluid particle acceleration 
in developed turbulent flow, we have shown that when the 
cross correlation parameter is taken in the exponential 
form, A ~ e' u ', is in a good qualitative agreement with 
the observed behavior of the experimental mean accel- 
eration conditional on velocity fluctuations u, as shown 
in Fig. This fit is of purely illustrative character per- 
formed with the only purpose to investigate effects pro- 
duced by nonzero cross correlation parameter A (it). We 
stress that the observed conditional mean acceleration 
is related to the flow anisotropy rather than to the cor- 
relation between stretching and vorticity controlled by 
A. Our primary purpose was to quantify to which extent 



variation of A affects the conditional variance of accelera- 
tion (a 2 \u). This effect has been found ignorably small as 
the result of the increase of additive noise intensity a(u) 
despite the fact that |A(it)| increases by about two orders 
of magnitude. We encountered ill-defined character of 
the distribution for big values of |u| that indicates either 
failure of the used numerical procedure or inappropriate- 
ness of this illustrative fit. For homogeneous isotropic 
case the conditional mean acceleration should be zero in 
contrast to the experimental data shown in Fig. For 
A = 0, the conditional acceleration distribution is well 
defined for any value of u. 

(ii) The additive noise intensity was taken to be a ~ 
e' u l that together with A ~ e' u have implied the varia- 
tion of the shape of the conditional probability distribu- 
tion function P{a\a{u), A(it)) which qualitatively agrees 
with the shapes of the experimental P(a\u) of the trans- 
verse component of acceleration at various values of u, as 
shown in Fig. ^ The increase of a tends to symmetrize 
the acceleration distribution. 

(iii) Variation of \u\ beyond certain value was found 
to imply ill-defined conditional distribution P(a\u) — 
P(a\a(u) , X(u)) probably because |A(u)| becomes compa- 
rable to the additive and multiplicative noise intensities 
which case requires a more detailed study of the nontriv- 
ial dependency of the distribution J2J on the parameters. 
If this is not due to a failure of the used numerical proce- 
dure at big | A | in intermediate calculations, this problem 
can be cured by the assumption that some of the other 
parameters of the model depend on u in certain way as 
well, to keep P(a\u) well-defined for any u. 
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